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Abstract: Gauge problem of monopole dynamics is studied in SU (2) lattice gauge theory. 
We study first abelian and monopole contributions to the static potential in four smooth 
gauges, i.e., Laplacian Abelian (LA), Maximally Abelian Wilson Loop (MAWL) and L- 
type gauges in comparison with Maximally Abelian (MA) gauge. They all reproduce the 
string tension in good agreement with the SU (2) string tension. MA gauge is not the only 
choice of the good gauge which is suitable for the color confinement mechanism. Using 
an inverse Monte-Carlo method and the blockspin transformation, we determine effective 
monopole actions and the renormalization group (RG) flows of its coupling constants in 
various abelian projection schemes. Every RG flow looks to converge to a unique curve 
which suggests gauge independence in the infrared region. 
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1. Introduction 

It is important to understand color confinement mechanism in Quantum Chromodynamics 
(QCD). Many numerical simulations have been done and they support the dual super- 
conductor scenario of the QCD vacuum as a confinement mechanism 0. Magnetic 
monopoles are induced by performing an abelian projection ||, i.e., a partial gauge-fixing 
that keeps 11(1)011(1). It is known that the string tension calculated from the abelian and 
the monopole parts reproduces well the original one when we perform an abelian projection 
in Maximally Abelian (MA) gauge where link variables are abelianized as much as possible. 
In addition to the string tension, many low energy physical properties of QCD are repro- 
duced from the abelian and the monopole degrees of freedom alone. It is called as "abelian 
and monopole dominance". These facts suggest that monopoles play an important role 
for the confinement mechanism. Actually, a low energy effective theory which is described 
in terms of monopole currents has been derived by Shiba and Suzuki Q and an almost 
perfect monopole action showing the scaling behavior is derived by Chernodub et al ||. 
Monopole condensation occurs due to energy-entropy balance Q. Abelian color-electric 
flux is squeezed into a string-like shape [||, by the superconducting monopole current. 
This squeezed color flux causes a confinement potential between quarks. 
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We note that we have infinite degrees of freedom when we perform an abelian pro- 
jection. That is to say, which gauge should be chosen? Recently Laplacian Abelian (LA) 
gauge was proposed and it appears to have similar good properties |8|, ||. Actually MA 
and LA gauges are very similar. Are MA and LA gauges exceptional? If such is the case, 
there must exist any reason to justify it, although it seems very difficult to find such a 
reason. Another interpretation is that monopole dynamics does not depend on the choice 
of gauge in the continuum limit, although it seems dependent of the gauge choice at the 
present stage of lattice study. In other words, MA gauge and LA gauge are considered to 
have a wider window even at present to see the continuum limit than other gauges have. 

Our aim of this paper is to show first that MA gauge is not a special choice of the 
good gauge for color confinement. We restrict ourselves to pure SU(2) QCD for simplicity. 
Here we discuss two new gauges in addition to LA gauge. They have a different continuum 
limit but they all can reproduce well the SU(2) string tension. The second is to derive 
an effective monopole action and to study the blockspin transformation of the monopole 
currents in various abelian projections. If their renormalization group (RG) flows converge 
onto the same line with a finite number of blockspin transformations, we can expect gauge 
independence of monopole dynamics in the infrared region. The paper is organized as 
follows: In Section 2, we present some theoretical and phenomenological arguments which 
support gauge independence of abelian and monopole dominance. In Section 3, we describe 
gauge fixing procedures being used. In Section 4, we show that the SU{2) string tension 
is well-reproduced from abelian or monopole degrees of freedom alone in four different 
abelian projection schemes. In section 5, we present our results from RG flow study of 
effective monopole actions in various abelian projections. In Section 6, we summarize our 
conclusions. 



2. Theoretical and phenomenological background 
2.1 Gauge fixings and abelian dominance 

It is known that the abelian Wilson loop reproduces well the SU(2) string tension numer- 
ically, if MA gauge or LA gauge is applied || |l(|. In the case of Polyakov gauge, the 
string tension which is calculated from abelian Polyakov loop correlators is exactly the 
same as that of SU(2) Shoji et al. have developed a stochastic gauge fixing method 

which interpolates between the MA gauge and no gauge fixing They have found that 
abelian dominance for the heavy quark potential is realized even in a gauge which is far 
from the MA gauge. In a finite temperature system, abelian Polyakov loops in various 
gauges reproduce the phase transition behavior of SU(2) Polyakov loop |13] ]. See, Figure 

B 

Abelian dominance is shown also analytically. Abelian Wilson loops constructed with- 
out any gauge fixing give the same string tension as that of SU(2) Wilson loops in the 
strong coupling expansion [10|. The same fact for any coupling region has been proved by 



-2- 



Ogilvie using the character expansion [14]. An abelian Wilson loop operator is given by 



W A [C] = -Tr 



n u v( s ) 



where is an abelian projected U(l) link variable. Since Wa is not a class function of 
SU(2) group, only the SU(2) invariant part extracted from Wa is non-vanishing in the 
expectation value. This can be written as 



wt = \\ v«\\ 



]J 9(s)u fl (s)g ] (s + fi) 

s./xGC 



Using a character expansion, we get an expression for the expectation value of the abelian 
Wilson loop in terms of SU (2) Wilson loops: 



P(C) 



{Wsu(2))l/2 + (half integer higher rep.). 



Since the lowest representation is dominant, we can show that the SU(2) string tension 
a SU(2) can be reproduced perfectly from the abelian string tension a a'- 

y , (W A (I+1,J+1))(W A (I,J)) 

aA = -i l ^J n (w A (i + i,J))(w A (i,J + i)) = asu ^ 

Furthermore, Ogilvie has shown that similar arguments hold even with a gauge fixing 
function 

if the gauge parameter A is small enough. 
2.2 Monopole dominance 

There are numerical results supporting monopole dominance. SU (2) string tension is well 
reproduced only from the monopole part of abelian Wilson loops in MA gauge |L5], and 
LA gauge ||. We note also that monopole Polyakov loops in various gauges reproduce the 
phase transition behavior of SU(2) Polyakov loop jl^l - See, Figure |2[ 

In addition to these numerical evidence, we can prove analytically gauge independence 



of monopole dominance if abelian dominance is gauge independent [17]. If abelian dom- 
inance is gauge independent, a common abelian effective action S e ff written in terms of 
abelian gauge field surely exists in any gauge and works well in the infrared region as in MA 
gauge. Since S e jf takes a form of modified compact QED, an effective monopole action 
can be derived analytically. One can evaluate the contribution of monopoles to the abelian 
Wilson loop using this effective monopole action. 

In MA gauge, it is known numerically that an effective monopole action composed 
of two-point self+Coulomb+nearest neighbor interactions is a good approximation in the 
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infrared region. The action can be transformed exactly into a modified compact QED 
action in the generic Villain form: 



Z = Pfl^exp - — (d6 + 2irn,AD(d6 + 2irn)) + i(J,0) 



where D ~ /3A _1 + a + 7 A. An expectation value of the abelian Wilson loop W = e 1 ^'^ 
can be estimated using this action, where J is the color electric current which takes a value 
±1 on a closed loop. When we use the BKT transformation [18, we get an expectation 
value of the abelian Wilson loop in terms of monopole currents k: 



(W) 



~ exp [-(k,Dk) -2ni(k,6A- 1 M) - ir 2 {J,(A 2 D)- l J) 



kez,dk=o 



where M takes a value ±1 on a surface whose boundary is J (J = 5M). Electric-electric 
current {J- J) interactions are of a modified Coulomb interaction and have no line singu- 
larity leading to a linear potential. A linear potential of the abelian Wilson loop originates 
from the second term of the monopole contribution. Gauge independence of monopole 
dominance is derived from that of abelian dominance. Gauge independence of an order 



parameter is also observed in Ref. [20|. 
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Figure 1: SU(2) Polyakov loop vs. abelian Polyakov loop in various gauges. The behavior of the 
SU(2) Polyakov loop is well reproduced by the abelian Polyakov loop in various gauges. 



2.3 An objection to gauge independence 

As we have shown in previous subsections, there is encouraging evidence which supports 
gauge independence of the confinement scenario in terms of monopoles. On the other hand, 
there is a strong objection to the idea of gauge independence. 

Consider a gauge called Polyakov gauge where Polyakov loop operators are diagonalized 



in the continuum finite-temperature QCD. It is proved [21, ^] that the singularities of the 



gauge fixing run only in the time-like direction. This means that there are only time-like 
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Figure 2: SU(2) Polyakov loop vs. monopole Polyakov loop in various gauges. The behavior of 
the SU(2) Polyakov loop is well reproduced by the monopole Polyakov loop in various gauges. 



monopoles in the system when the Polyakov gauge is employed, if the degeneration points 
in abelian projection only correspond to monopoles as 't Hooft argued. Since such time-like 
monopoles do not contribute to the physical string tension p3[| , monopole dominance is 
violated. 

But numerically the above theoretical expectation seems to be inconsistent with nu- 
merical data. We show our preliminary result in Figure ||[ The spatial and temporal 
monopole densities are plotted in Figure | as a function of lattice spacing a in the unit of 
physical string tension y/dp. These densities are defined as 

n (R\- l^sEi=l,2,3\ k i( S )\ „ fg\ _ E,|fe*00l 

respectively. Figure |3| shows that spatial (lattice) monopole density may take non-zero 
value even in the a — > limit. This is not compatible with the theoretical expectation 
above. In the authors' opinion, the continuum limit of lattice monopoles must contain 
extra ingredients different from the expected monopoles corresponding to singularities of 
Polyakov loop operators. We will give a detailed analysis elsewhere. 

3. Various abelian projections on a lattice 

To check gauge (in) dependence of monopole dynamics, we study the abelian projection in 
various gauges. 

1. MA gauge: 

The most known is Maximally Abelian (MA) gauge. It is defined by maximizing the 
following quantity Rma- 



Rm a = Tr Y,M s )^Ul(s)a 3 . (3.1) 



8,/J. 
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Figure 3: Monopole density in Polyakov gauge versus lattice spacing 

This is achieved by diagonalizing an operator 

Xma(s) = (yMoJjfa) + Ul(s - ftoMs - A)) • 



That is, 



Xma(s) -> X' MA {s) = V(s)X M a(s)V\s) = diag{Ai,A 2 }, 



where V(s) is a gauge transformation matrix. The diagonalization corresponds to 
the condition, 

J2{d»TiAl)A± = 0, (3.2) 

in the continuum limit. 
2. LA gauge g: 

First consider MA gauge again. To maximize Rma hi Eq.(|3.1|) is to minimize the 
functional 



(3.3) 



Sma = ^{l-^TV[$( S )C/ /1 ( S )$(5 + A)^J(s)] 

s,fj, *■ 

= ^{^(,)<( S )/( S + A)}, 

S,fJ, 

where is a gauge field in the adjoint representation, 

<(«) = ^Tr (cTaU^UKs)) . 
is parameterized by a spin variable <p which satisfies a local constraint 

3 3 

^(s) = vHs)a 3 V(s) = Y j nsW , £(« 2 = 1. 



(3.4) 



a=l 



0=1 
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Because of the local constraint from the normalization, it is very difficult to find a 
set of (j) which realizes the absolute minimum of Eq.(3.3). 

The key idea of the LA gauge fixing is to relax this constraint: 

£fo») 2 = i — EX> a (*)) 2 = 1 

a=l s a=l 

The functional to minimize becomes 



where 



2 



~ a ty = E ( 2 V - <(*)<W - - A)<WJ • (3-6) 

Minimizing Eq. ( |3.5| ) amounts to finding the eigenvector belonging to the lowest eigen- 
value of the covariant laplacian operator. This eigenvalue problem can be solved 
numerically (we used an implicitly restarted Arnoldi method. For example, see 
Ref. |24||). The gauge transformation matrix V(s) is defined by 

3 

vHs)a 3 V(s) = J24> a (s)va, (3.7) 
o=l 

where 

3 

ns) = P (sW(s) , P 2 ( S ) = Y,(r(s)) 2 . (3.8) 

a=l 

In the continuum limit, LA gauge corresponds to the gauge condition, 

(d, T iA%) (p 2 A±) = 0. (3.9) 



3. MAWL gauge J2§: 

Maximally Abelian Wilson Loop (MAWL) gauge is a gauge which maximizes a Wilson 
loop operator written in terms of abelian link variables: 

W A = cos 9^(8), (3.10) 

where 0^(s) = O^s) + 9 v (s + p) — 9^(s + v) — 9 u (s). It is achieved by diagonalizing 
the following operator 

sin^(s) - sin 0^(5 - it) 



Xmawl(s) = E 



U 2 (s,fi) + Ui( S ^) 



-(U(s,fi)a 3 U\s,fi)) 



smd^s-n-u )-sm0 lu , (s - /i) t A A 

+ pT^T : — x , TT 2( — — \ ( u (s-h,h)<t 3 U (s-fi,n ) 

Utf(s- n,fj,) + U${s- fj,,fx) 

In the continuum limit, we get the following gauge condition: 

J2duU„A± = 0. (3.11) 
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4. L-type gauge: 

There are infinitely many gauges similar to MA gauge. Here we show one simplest 
extension called L-type gauge. It is defined by maximizing 

Rl = Tr U^{s)a 3 U u {s + fi)a 3 Ut(s + p)a 3 Ut(s)a 3 . (3.12) 
This is given by diagonalizing 

+Ul(s - fi)a 3 U v {s - iX)a 3 Ul{s - (i)a 3 U^{s - fi) 

A schematic representation of Rl is shown in Figure |j. 

In the continuum limit, we get the following gauge condition: 

]T{(<9 M ± iagAl) + (d v ± iagAl)}{Aj + A$) = 0. (3.13) 

5. There are various gauges called unitary gauge. Polyakov gauge and Fj.2 gauge are 
defined with the following operators which are diagonalized: 

xpoi(s) = n u ^ s + (* - !) a )' ( 3 - 14 ) 

i=l 

X Fl2 (s) = Ui(s)U 2 (8 + l)ul(s + 2)t4(s), (3.15) 

respectively. 

In the continuum, the Polyakov gauge is reduced to 

A±(x) = 0, (3.16) 

whereas F12 gauge gives 

F 1 %(x) = 0. (3.17) 

6. We also consider simple abelian components extracted without gauge-fixing, where 



exact abelian dominance is proved analytically [14]. 



4. String Tension 

As a first step, we measure abelian and monopole contributions to the string tension 
in various abelian projections. We used 100 configurations of 32 3 x 16 lattice for the 
measurement. In this case, the critical point lies near /3 C ~ 2.7. We set the gauge coupling 
(3 to 2.5, so that the system is in the confinement phase. To reduce the statistical errors 



efficiently, we adopted the hypercubic blocking [26] to the original configurations. 
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The value of Polyakov loop correlators corresponds to the 
static potential between one pair of quark and anti-quark; 

(TrP(0)TrP t (.R)) = e~ viR)/T , (4.1) 

where P(R) is the Polyakov loop operator Eq.( |3.14] ). V(R) gives 
the inter-quark potential 

V(R) = aR-^-+c 
R 

and T = l^N^a) is the temperature of the system. 
The abelian Polyakov loop operator is written as 

JV4-I 

P a = exp[i ^ A {s + il)}. 



(4-2) G3 



a 



Figure 4: Schematic 
representation of L-type 
gauge 

(4.3) 



Eq.([4,3j) can be decomposed to photon and monopole parts [11] as follows: 



Pa 



P ■ P 

1 p 1 mi 



2V4-1 



P p = expH Y, D ^ +il - s '^^ s '^ 

i=0 s> 
N 4 -l 

P m = exp[-27n V J" D(s + il - s')d' u n uA (s')}, 



i=0 s' 



where we use an identity 



s> 



D(s-s') [d' u @ vi {s') + d^d'Ms'))] . 



The abelian field strength tensor 

e^(s) = ^(s)+^(s + /i)-^(s + z>)-^(s) , (-Air < M s ) < 47r ) 

can be decomposed into two parts: 

Qnu(s) = Qfj.uis) + 27rn^(s) , (-ir < 6^(s) < ir). 

Here, Q^ u (s) is interpreted as the electro-magnetic flux through the plaquette and the 
integer valued n^ v (s) corresponds to the number of Dirac string piercing the plaquette. 
D(s — s') is the Coulomb propagator on a lattice. 

Figures ||, ||, [7| and || show the values of SU(2), abelian and monopole Polyakov loop 
correlators in MA, LA, MAWL and L-type gauges, respectively. The values of abelian and 
monopole Polyakov loop correlators in each gauge almost degenerate. The string tension a 
can be extracted from these values by fitting them to Eq.( |4.1| ). Fitted lines are also plotted 
in the same figure. In the case of MA gauge, fitted values are consistent with the results 
by Bali et al p7[ |. In the case of other gauges like a unitary gauge, one can not extract the 
string tension clearly from the abelian and the monopole Polyakov loop correlators due to 
large statistical errors. 

Explicit values of the fitted string tension are shown in Table |l| They almost agree each 
other, although these four gauges have different gauge fixing condition in the continuum 
limit. 
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MA 


LA 


MAWL 


L-type 


Abelian 


0.03054(45) 


0.03011(34) 


0.03051(45) 


0.03065(43) 


Monopole 


0.02545(31) 


0.02536(28) 


0.02546(31) 


0.02624(34) 



Table 1: Fitted string tensions (32 3 x 16 lattice, (3 = 2.5) <7 SU ( 2 ) = 0.03446(105). 



5. RG flows of the effective action in various abelian projections 

To clarify what is happening in the monopole dynamics, we study the effective monopole 
actions in various gauges in this section. 



5.1 Simulation method 

Our method to derive an effective monopole action is the following. We generate SU (2) 
gauge fields {U^(s)} using the standard SU(2) Wilson action. We consider 48 4 hypercubic 
lattice for [3 from 2.1 to 2.5. We took 50 independent configurations after 10,000 thermal- 
ization sweeps. Then, we perform the abelian projection in six different gauge fixings to 
extract abelian gauge fields {u^{s)} from SU{2) gauge fields. 

One can define magnetic monopole currents from abelian field strength tensor along 



* c — i 1 1 1 1 1 1 , 1 




Distance R Distance R 

Figure 5: Abelian and monopole Polyakov Figure 6: Abelian and monopole Polyakov 
loop correlator in MA gauge loop correlator in LA gauge 
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the way of DeGrand and Toussaint [28]. We can define the monopole current k tJi (s) as 



M S ) = ^iivpadyTlp^S + £). (5.1) 



Vu 



By definition, it satisfies a current conservation law 

d'M s ) = °> 

where and 3^ denote the forward and the backward differences in /i-direction respec- 
tively. 

We want to get an effective monopole action S[k] on the dual lattice, integrating out 
all degrees of freedom except for monopoles: 

Z = j VUe-^^iX^ApiU) 

Vue~ s eff [u] 
JIT.) J Vue~ s <ffM5(k,f(u)) 
= 11 E (U^M,o)e-^, 

fe M (s)=-oo \m,v / 

where Un = C^u^ and X 1 * 1 is the off-diagonal element of the matrix X which is diagonalized 
in the procedure of abelian projection. Ap(U) is the Faddeev-Popov determinant and 
8(k,f(u)) gives the definition of the monopole current k as a function of abelian gauge 
field u. 

Above integrations are done numerically. We create vacuum ensembles of monopole 
currents using the Monte-Carlo method and the definition of the monopole current Eq.(|5.l|), 
Then, we construct the effective monopole action from monopole vacua using the Swend- 



sen's inverse Monte-Carlo method which was developed originally by Swendsen [29] and 
extended by Shiba and Suzuki ||]. 

We consider a set of independent and local monopole interactions which are summed 
up over the whole lattice. We denote each interaction term as Si[k]. Then the effective 
monopole action can be written as a linear combination of these operators: 

S[fc]= J>Si[fc], (5.2) 

i 

where gi denotes the effective coupling constants. Explicit forms of the interaction terms 
are listed in Table § and ||. We determine the set of couplings from the monopole 
current ensemble {k^(s)} with the aid of an inverse Monte-Carlo method. Practically, we 
have to restrict the number of interaction terms. The form of action adopted here is 27 
quadratic interactions and 4-point and 6-point interactions |B|, 30]. 
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We perform a blockspin transformation in terms of the monopole currents on the dual 
lattice to study the RG flow. The n-step blocked current is defined by 

n-1 

K M n) )= Yl ^(ns^ + (n-l)p + iO + jp + la). (5.3) 

i,j,l=0 

The blocked lattice spacing b is given as b = na{(3) and the continuum limit is taken as the 
limit n — > oo for a fixed physical scale b. We determine the effective monopole action from 
the blocked monopole current ensemble {K fl (s^)}. Then one can obtain the RG flow in 
the 29-dimensional coupling constant space. 



Coupling 


Distance 


Type 


Coupling 


Distance 


Type 


9i 


(0,0,0,0) 




915 


(2,1,1,0) 


kn(s + 2p + + p) 


92 


(1,0,0,0) 


kfi(s + p) 


9m 


(1,2,1,0) 


k^(s + p + 20 + p) 


93 


(0,1,0,0) 


k^(s + 0) 


917 


(0,2,1,1) 


k^s + 20 + p + a) 


94 


(1,1,0,0) 


k^(s + p + 0) 


918 


(2,1,1,1) 


k^s + 2p + + p + a) 


55 


(0,1,1,0) 


k^(s + + p) 


919 


(1,2,1,1) 


k^s + p + 20 + p + a) 


96 


(2,0,0,0) 


kfi(s + 2p) 


920 


(2,2,0,0) 


k^(s + 2p + 20) 


97 


(0,2,0,0) 


k^(s + 20) 


921 


(0,2,2,0) 


k^(s + 2v + 2p) 


9s 


(1,1,1,1) 


kfj,(s + p + + p + a) 


922 


(3,0,0,0) 


k^(s + 3p) 


99 


(1,1,1,0) 


M s + p + + p) 


923 


(0,3,0,0) 


k^s + 30) 


9w 


(0,1,1,1) 


kfj,(s + + p + a) 


924 


(2,2,1,0) 


k^s + 2p + 20 + p) 


9n 


(2,1,0,0) 


k^{s + 2p + 0) 


925 


(1,2,2,0) 


k lt (s + fi + 2u + 2p) 


912 


(1,2,0,0) 


k^{s + p + 20) 


926 


(0,2,2,1) 


k^s + 20 + 2p + a) 


913 


(0,2,1,0) 


k^s + 20 + p) 


927 


(2,2,1,0) 


k p (s + 2p + 20 + p) 


9u 


(2,1,0,0) 


k u (s + 2p, + 0) 









Table 2: The quadratic interactions used for the modified Swendsen's method. 



5.2 Numerical Results ^ v m 

Coupling lype 



The effective monopole action is determined A „ •„+ „ /v^4 i,2/„\ 

4-pomt g 2 8 z^s { 2^^=-4>i s i 

successfully. All coupling constants which ). 

are contained in the effective monopole ac- 

6-point 929 E s ( E M =-4 k U s ) 

tion are obtained with relatively small er- Table 3 . The higher Qrder interactions used for 
rors. We use the jackknife method for the thc modified Swendsen's method, 
error estimation. These effective monopole 

actions except in MA gauge are determined for the first time in this paper. Moreover, these 
effective monopole actions are determined from the blocked monopole configurations, too. 
The results are summarized as follows: 

1. Only the quadratic interaction subspace seems sufficient in the coupling space for the 
low-energy region of QCD. Figure |9] and 1C show coupling constants 1 for 4-point and 



1 Effective coupling constants for the blocking factor n = 1 are omitted in Figures ^| and |l4|-| 
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6-point interaction terms versus physical scale b. In the case of MA, LA, MAWL and 
L-type gauges, these coupling constants take relatively larger absolute values for small 
b region. They become negligibly small for large b region. In the case of Polyakov, 
F12, no gauge fixings, coupling constants for 4-point and 6-point interaction terms 
take the values very close to zero in the whole region of b. 



2. Typical case of the coupling constants for quadratic interaction terms versus squared 



distances in the lattice unit are shown in Figure 11. We see that coupling constants 
for the self interaction term g\ and the nearest-neighbor interactions 52 and 53 are 
dominant, and 52 — 93- Other couplings decrease exponentially as distance between 
the two monopole currents grows. Such a behavior does not depend on a gauge 
coupling constant j3. Therefore, we concentrate our analysis on the coupling constants 
of quadratic interaction terms, especially g\ and g<i- 



3. We used a standard iterative gauge fixing procedure for MA, MAWL and L-type 
gauges. In such a case, gauge fixing sweeps may be stuck for some local minima of a 
gauge fixing functional. Different local minima give rise to different gauge transfor- 
mations, but they can not be distinguished from the view point of the iterative gauge 
fixing procedure. These are the lattice Gribov copies. Indeed, Bali et al. showed that 
the effect of such copies to the abelian string tension is not so small |27| . To check the 
effect of copies to the effective couplings, we generate 100 of SU(2) configurations on 
24 4 lattice at (3 = 2.5. Then, we generate 7 of gauge equivalent configurations (i.e., 
copies) via a random gauge transformation. Using these gauge copies, we constructed 



effective monopole actions and compared their effective couplings. Figure 12 shows 
gi in the case of MA gauge. g\ for the different blocking factors are described in 
different symbols. We see some fluctuations in g\ for MA gauge. This is nothing but 
the effect of lattice Gribov copies. The effect of the copies, however, are negligibly 
small. Therefore, qualitative analyses which are given later will not be affected. In 
principle, LA gauge does not have such a copy ||]. Indeed, we confirmed that effective 
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Figure 11: Effective couplings vs. squared distances in lattice unit. (MA gauge, (3 — 2.1,2.3 and 
2.5, effective couplings for n — 8 blocked monopole) 



couplings for LA gauge are not affected by Gribov copies (Figure |l3|) 



4. Figure 14 and Figure 15 show the most dominant quadratic self coupling constant 
gi and quadratic nearest-neighbor coupling constant g2 versus physical scale b in the 
case of MA, LA, MAWL and L-type gauges, respectively. In these gauges, effective 
coupling constants take large values in small b region and the scaling behavior (i.e., 
a unique curve for different blocking factor n) is seen even in small b region. The 
effective actions which are obtained here appear to be a good approximation of the 
action on the renormalized trajectory corresponding to the continuum limit. In ad- 
dition to this, coupling constants for these four gauges are very close to each other, 
although these gauges have a completely different form in the continuum limit. 
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Figure 12: Gribov copy effect for g\ (MA 
gauge) 
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However, in the case of Polyakov, F12 and no gauge fixings, coupling constants are 
different from those in the above four gauges (See, Figure |0] and Figure |l7| ). In 
these gauges, coupling constants take smaller values and the scaling behavior is not 
seen especially in small b region. To clarify the scaling properties of these coupling 
constants, we show the figures as showing a distinction between the different blocking 
factors n in two typical gauges. In the case of Polyakov gauge (Figure 18), the 
coupling constants depend on the blocking factor n strongly in small b regions. On 
the other hand, in the case of LA gauge (Figure IS), renormalized coupling constants 
are lying on the unique curve. 



6. Once the effective actions are fixed, we can see from energy-entropy balance of the 
system whether monopole condensation occurs or not. If the entropy of a monopole 
loop exceeds the energy, the monopole loop condenses in the QCD vacuum. In 
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Figure 14: The most dominant self coupling 
gi vs. physical scale b in MA, LA, MAWL 
and L-type gauges. 



Figure 15: Nearest-neighbor coupling 32 vs. 
physical b in MA, LA, MAWL and L-type 
gauges. 
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Figure 16: The most dominant self coupling 
gi vs. physical scale b in MA, Polyakov, F12 
and no gauge fixings. 



Figure 17: Nearest-neighbor coupling 52 vs. 
physical scale b in MA, Polyakov, F12 and no 
gauge fixings. 



four-dimensional lattice theory, the entropy of a monopole loop can be estimated as 
In 7 per unit loop length. It is determined by the random walk without backward 
tracking. The action can be approximated by the self interaction term g\ alone since 
the interactions with two separate currents are almost canceled [31]. The free energy 
per unit monopole length is approximated by 



F 



91 



In 7, 



since g\ can be regarded as the self energy per unit monopole loop length. If g\ < In 7, 
the entropy dominates over the energy, that is, monopole condensation occurs. In 
Figure 14 and Figure we see that the entropy of the system dominates over the 
energy in the large b region for all gauges. In other words, monopole condensation 
occurs in the large b region for all gauges. 



Figures 2C ,21,22 and 23 show the RG flows projected onto gi-g2, ffi-<?5> 9\-Qi and 
gi-gio coupling planes, respectively. The effective coupling constants for all gauges 
seem to converge to the identical line for large b region. This may show gauge 
independence of the monopole condensation in the low energy region. Although all 
coupling constants become very small in the large b region, it is important that the 
slopes of the renormalization flows seem to converge in all gauges. 



6. Summary 

We have measured first the abelian and the monopole contributions to the string tension 
in four types of abelian projections, i.e., MA, LA, MAWL and L-type gauges. They show 
a good agreement with each other. Monopole string tension are extracted in the same 
manner as abelian string tension, and they agree also with each other. MA and LA gauges 
are not unique good gauges. 

Next, we have determined the effective monopole actions in various gauges from 
monopole vacua using the modified Swendsen's method. In the case of MA gauge, an 
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Figure 18: g\ versus b in MA gauge and 
Polyakov gauge. Each symbols correspond 
to the different blocking factors n. 
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effective monopole action has already been obtained in Ref. In addition to this action, 
the effective monopole actions in Polyakov gauge, F12 gauge, LA gauge, MAWL gauge, 
L-type gauge and no gauge fixing are also determined for the first time in this paper. 
Moreover, these effective actions are determined on the blocked monopole vacua, too. In 
these effective actions, two point interactions are dominant, whereas 4-point and 6-point 
effective coupling constants are negligibly small in the infrared region. The RG flows seem 
to converge to the identical line when repeating the blockspin transformation. It is im- 
portant that the slopes of renormalization flows in all gauges seem to converge. The data 
are compatible with the assumption of gauge independence of monopole dynamics in the 
continuum limit. Energy-entropy balance also tells us the monopole condensation occurs 
in the large b region for all gauges. 
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Figure 20: RG flows of the effective 
monopole actions project onto the 51-172 cou- 
pling plane 



Figure 21: RG flows of the effective 
monopole actions project onto the 51-35 cou- 
pling plane 
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A. Maximally Abelian Wilson Loop (MAWL) gauge 

SU(2) gauge field U^(s) can be parameterized by its isospin components. In this section, 
we denote each isospin component of U^(s) as Uo(s,fi), Ui(s,fi) and so on, for simplicity. 
This gauge is realized with maximizing the abelian Wilson loop of 1 x 1 size: 

R= Yl cos0 ^( s )> (A.l) 

where the abelian link field is extracted as 

0(s, n) = arctan([/ 3 (s, fi)/U (s, fi)). (A.2) 
Let us consider an infinitesimal gauge transformation of U, 

U'(s,n) = (1 + iai(s)ai)(U (s,fj,)I + iUj(s,fj,)(Tj)(l - ia k (s + fi)a k ). (A.3) 

This gives 

5U (s, /i) = -(oi(s) - ai(s + fi))Ui(s, n), (A.4) 
5U k (s,n) = (ak{s)-a k {s + p))Uo(s,fx)-eijk(ai(s) + ai(s + p,))Uj(s,fi). (A.5) 

Then R changes as 

SR = - Y snQ^(s)(S9(8,ii) +66(s + fl,v) -S9(s + fl,fi) -60(s,v)), (A.6) 



dV{s,fi) = 2 2 . (A.7J 



where 

U (s, fi) SU 3 (s,fx) - U 3 (s,fi)5 U (s, fi) 

One can check that R is invariant under the U (1) transformation. Hence we do not 
need to consider the 03(5) part. First, let us consider the a\ part. Since there is the sum 
over whole lattice sites s, one can shift the site variable, for example, s to s — ft. Also one 
can use the (anti) symmetric property with respect to fi and v directions. Finally one gets, 

— 2~ = Yl ( a i( s ) x i( s ^^) +a 2 (s)X 2 (s,fi,is)), 

Y( x U x {8,n)Uz{8,n) - U (s,n)U 2 (s,n) 
Xi(s,n,v) = Wi 2 2 — — — — 



+w 2 



U*(s,») + Ui(s,ti) 
Ui(s - A, n)U${s - A, At) + Up(s - /t, n)U 2 (s - £1, n) 
Uq(s- fi,ii) + Ul(s- fan) 



- 18 - 



v , B , ,a U / U 2 (s,fj l )U 3 (s,fj 1 ) + Uo{s,n)U 1 (s,fj,) 

A 2 (S, /X, i/J = Wi 



+ W 2 



£/ 2 (s,M) + f7 3 2 (s,/x) 

J7 2 (S - p,,/J,)U 3 (s -fl,fJL)- U (s - p,,fjL)Ul(s - ft,/ji) 



U$(s- fi,iJ,) + Ui(s- fafj,) 
Wi = sinS^fs) - sin6 M „(s - z>), 

W 2 = sin0 MI/ (s - /t- z>) - sin0 Mi ,(s - /t). 

When we write X 1 * 1 = Xi ± iX 2 , it is easy to see X 1 * 1 transforms covariantly under the 
residual U(l). 

Finally, one gets the matrix which is diagonalized in this gauge, 



sm 8^(s) - sin 6^Qs - t/) , 
U (s,fi )+Ui f(s,/i) 



+ 7727 : — T— 7727 — x (s - fi, fi)a 3 U (s-/i,/i ) 

Because of the non-locality of the gauge condition, one can not calculate the gauge trans- 
formation matrix which diagonalizes X(s) in a simple way. Therefore, we employed an 
iterative updation procedure to satisfy the gauge condition. 

1. Make a trial gauge transformation, adopting a>\ and «2 as follows: 
ol\(s) = —kXi(s), 

0:2(5) = -kX 2 (s). 

2. Measure R. If R becomes larger than before, accept this trial and repeat step 1. If R 
becomes smaller than before, take K new = K id/2 and adopt the gauge transformation 
using this K new with respect to the configuration before trial, and then repeat step 1. 

3. If the off-diagonal element of X(s) becomes smaller than a suitable threshold (we set 
this to 1.0), one can regard the gauge fixing procedure as having been completed. 

We set an initial value of k to 0.1. R can be maximized as long as we take k > 0. 
We apply the MA gauge fixing as a preconditioning for MAWL gauge fixing and then 
we perform the above procedure to the MA fixed configuration. This preconditioning is 
required to improve a convergence property of the MAWL gauge fixing. We have to note 
that the configurations which obtained via the above procedure are not perfectly gauge 
fixed because the off-diagonal element of X(s) still remain not very small. 
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